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Optimized Uplink Transmission in Multi-Antenna 
C-RAN with Spatial Compression and Forward 

Liang Liu and Rui Zhang 

Abstract 

Massive multiple-input multiple-output (MIMO) and cloud radio access network (C-RAN) are two promising 
techniques for implementing future wireless communication systems, where a large number of antennas are deployed 
either being co-located at the base station (BS) or totally distributed at separate sites called remote radio heads 
(RRHs) to achieve enormous spectrum efficiency and energy efficiency gains. In this paper, we consider a general 
antenna deployment design for wireless networks, termed multi-antenna C-RAN, where a flexible number of 
antennas can be equipped at each RRH to more effectively balance the performance and fronthaul complexity 
trade-off beyond the conventional massive MIMO and single-antenna C-RAN. In order to coordinate and control 
the fronthaul traffic over multi-antenna RRHs, under the uplink communication setup, we propose a new “spatial- 
compression-and-forward (SCF)” scheme, where each RRH first performs a linear spatial filtering to denoise and 
maximally compress its received signals from multiple users to a reduced number of dimensions, then conducts 
uniform scalar quantization over each of the resulting dimensions in parallel, and finally sends the total quantized 
bits to the baseband unit (BBU) via a finite-rate fronthaul link for joint information decoding. Under this scheme, 
we maximize the minimum signal-to-interference-plus-noise ratio (SINR) of all users at the BBU by a joint resource 
allocation over the wireless transmission and fronthaul links. Specifically, each RRH determines its own spatial 
filtering solution in a distributed manner to reduce the signalling overhead with the BBU, while the BBU jointly 
optimizes the users’ transmit power, the RRHs’ fronthaul bits allocation, and the BBU’s receive beamforming with 
fixed spatial filters at individual RRHs. Through numerical results, it is shown that given a total number of antennas 
to be deployed, multi-antenna C-RAN with the proposed SCF and joint optimization significantly outperforms both 
massive MIMO and single-antenna C-RAN under practical fronthaul capacity constraints. 

Index Terms 

Multi-antenna cloud radio access network (C-RAN), massive multiple-input multiple-output (MIMO), fronthaul 
constraint, spatial-compression-and-forward (SCF), power control, beamforming, signal-to-interference-plus-noise 
ratio (SINR) maximization. 


I. Introduction 

It is anticipated that the mobile data traffic will grow 1000 times higher in the next decade with a rate of roughly 

a factor of two per year |[T). As a result, the fifth-generation or 5G wireless system on the roadmap has to make a 

major paradigm shift to accommodate this dramatic explosion in future demands for wireless data communications. 

Both as highly promising candidate techniques for the 5G wireless network, massive multiple-input multiple-output 
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Fig. 1. Comparison of massive MIMO, single-antenna C-RAN, and multi-antenna C-RAN in the uplink transmission. 


(MIMO) 0, 01 and cloud radio access network (C-RAN) 0] have recently drawn significant attentions (see Fig. 
[I]). Both techniques advocate the use of very large number of antennas and centralized signal processing to achieve 
enormously improved efficiency in spectrum and energy usage. However, they are practically implemented with 
different performance and complexity trade-off considerations. 

In a massive MIMO system, as depicted in Fig. Q2a), a large number of antennas are deployed at the base station 
(BS) to reap all the benefits of the conventional multi-user MIMO system but with a much greater scale. It is shown 
in j2l that in the asymptotic regime where the number of antennas at the BS is much larger than that of the users, 
the channels of different users are orthogonal with each other, and thus a simple matched filter, i.e., maximal ratio 
combining (MRC) in the uplink and maximal ratio transmission (MRT) in the downlink, is capacity-achieving. 
Moreover, the required transmit energy per bit vanishes as the number of antennas goes to infinity {5|. Despite of 
the above appealing benefits, there arc many issues that potentially limit the performance of massive MIMO. For 
instance, due to the limited space for installation, the number of antennas at the BS is finite and their channels 
are generally correlated; as a result, the theoretical performance limit in the asymptotic regime cannot be fully 
achieved in practice with simple MRC/MRT. Moreover, similar to conventional MIMO system, the performance of 
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the cell-edge users in the massive MIMO system is still a bottleneck since their channel conditions may be weak 
due to more significant path attenuation. 

On the other hand, in a single-antenna C-RAN as shown in Fig. [Ilb), each antenna is deployed at one separate 
site called remote radio head (RRH), which is connected to one central unit (CU) via a high-speed fronthaul link 
(fiber or wireless) for joint information processing. Unlike the BS in the massive MIMO system which decodes 
or encodes the user messages locally, in C-RAN each RRH merely forwards the signals to/from the CU via its 
fronthaul link, while leaving the joint encoding/decoding complexity to a baseband unit (BBU) in the CU. As 
a result, if the data can be perfectly transmitted between the BBU and RRHs over the fronthaul links without 
any error, the single-antenna C-RAN can be viewed as another form of massive MIMO but with the antennas 
densely distributed over the whole served region such that all the mobile users can be served by some adjacent 
antennas with strong channel conditions. However, one obstacle constraining the practically achievable throughput 
of C-RAN is the limited-capacity fronthaul link between each RRH and BBU. According to H, the capacity 
of the current commercial fibers, which is about several Gbps (gigabits per second), can be easily overwhelmed 
in C-RAN even under moderate data traffic; whereas wireless fronthaul techniques such as those exploiting the 
millimeter-wave bands in general operate with even lower data rate per link and shorter distance as compared to 
the fiber based solution. In the literature, a considerable amount of effort has been dedicated to study fronthaul 
quantization/compression techniques in the uplink communication ll6l- llT2l . Specifically, the “quantize-and-forward 
(QF)” scheme |[T3l - llT5l is adopted to reduce the communication rates between the BBU and RRHs, where each 
single-antenna RRH samples, quantizes and forwards its received signals to the BBU over its fronthaul link. 
Moreover, the QF scheme is usually studied under an information-theoretical Gaussian test channel model, and 
the quantization noise levels of all the RRHs are jointly optimized to maximize the end-to-end throughput subject 
to the capacity constraints of individual fronthaul links. Recently, llT6l reveals that subject to the capacity of the 
current commercial fibers, the above theoretical performance upper bound can be approached by applying a simple 
uniform scalar quantization at each RRH for orthogonal frequency division multiplexing (OFDM) based single¬ 
antenna C-RAN. Moreover, another line of research in C-RAN studies the control of the overall fronthaul traffic 
by selectively activating only a small subset of all the RRHs to meet the user demands ifTTl . 

In the aforementioned massive MIMO system and single-antenna C-RAN, the antennas are deployed in a totally 
co-located and separated manner, respectively. In addition to these two extreme cases, a more general solution is 
to distribute the antennas at a reduced number of sites as compared to single-antenna C-RAN, while the RRH 
at each site is equipped with multiple antennas. This so-called multi-antenna C-RAN design can provide both 
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improved channel conditions to the users as in single-antenna C-RAN as well as local spatial multiplexing gain 
as in massive MIMO; furthermore, by adjusting the number of antennas at each RRH and the number of RRHs, 
these benefits can be flexibly traded-off with a given total number of antennas deployed. However, the fronthaul 
issue becomes more severe in the multi-antenna C-RAN as compared to conventional single-antenna counterpart, 
since the limited fronthaul capacity at each RRH needs to be allocated over the outputs of more antennas, thus 
leading to more quantization errors in general. To tackle this challenge and address the above design trade-off, in 
this paper we focus our study on the uplink communication in C-RAN with multi-antenna RRHs each subject to 
a finite fronthaul capacity constraint. The contributions of this paper are summarized as follows. 

• We propose a novel “spatial-compression-and-forward (SCF)” scheme for the multi-antenna RRHs to balance 
between the information conveyed to the BBU and the data traffic over the fronthaul links. Specifically, 
each RRH first performs spatial filtering to its received signals at different antennas from the users such that 
the signals are denoised and maximally compressed into a reduced number of dimensions. Then, the RRH 
applies the simple uniform scalar quantization lfl6l over each of these dimensions in parallel by appropriately 
allocating the number of quantization bits over them under a total rate constraint. Last, the quantized bits by 
RRHs are sent to the BBU for joint information decoding via their fronthaul links. It is worth noting that 
at each RRH, the compression of the correlated signals received by all its antennas is fulfilled by a simple 
linear filter, rather than the complicated distributed source coding, e.g., Wyner-Ziv coding, as in single-antenna 
C-RAN where all the antennas are physically separated lUTl . liT2l . 

• With the proposed SCF scheme, we formulate the optimization problem of joint users’ power allocation, RRHs’ 
spatial filter design and quantization bits allocation, as well as BBU’s receive beamforming to maximize the 
minimum signal-to-interference-plus-noise ratio (SINR) of all the users in the uplink. We propose an efficient 
solution to this complicated design problem for ease of practical implementation. First, each RRH computes 
its own spatial filter in a distributed manner based on its locally received signal covariance matrix. This 
distributed solution helps save significant fronthaul resources for exchanging control signals with the BBU. 
Then, the BBU jointly optimizes the other parameters to maximize the minimum SINR of all users based 
on the alternating optimization technique. Specifically, with given quantization bits allocation at each RRH, 
we extend the well-known fixed-point method lfT8l to obtain the optimal transmit power levels for the users 
and receive beamforming vectors at the BBU. On the other hand, given the above optimized parameters, the 
corresponding quantization bits allocation at each RRH is efficiently solved. 

• Last, we investigate the following interesting question: Given a total number of antennas for a target area, 
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should all of them be deployed at one BS, i.e., massive MIMO, or one at each RRH, i.e., single-antenna C- 
RAN, or optimally divided over a certain number of RRHs, i.e., multi-antenna C-RAN, given a practical total 
fronthaul capacity constraint for C-RAN? Through numerical examples, it is shown that with the proposed SCF 
scheme and joint wireless-fronthaul-BBU optimization, multi-antenna C-RAN generally performs noticeably 
better than both massive MIMO and single-antenna C-RAN under moderate fronthaul rate constraints. 

The rest of this paper is organized as follows. Section HD presents the system model for C-RAN with the 
SCF scheme. Section Hill formulates the minimum SINR maximization problem. Section [TV] presents the proposed 
solution for this problem. Section [V] provides numerical results to verify the effectiveness of the proposed multi¬ 
antenna C-RAN with SCF. Finally, Section [VI] concludes the paper. 

Notation : Scalars are denoted by lower-case letters, vectors by bold-face lower-case letters, and matrices by bold¬ 
face upper-case letters. I and 0 denote an identity matrix and an all-zero matrix, respectively, with appropriate 
dimensions. For a square matrix S, S A 0 (S A 0) means that S is positive (negative) semi-definite. For a matrix 
M of arbitrary size, M u and rank (S') denote the conjugate transpose and rank of M, respectively, diag(a) 
denotes a diagonal matrix with the main diagonal given by the vector a, while diag(Ai, • • • , Ak) denotes a block 
diagonal matrix with the diagonal given by square matrices A\, ■ ■ ■ , Ak- E[f] denotes the statistical expectation. 
The distribution of a circularly symmetric complex Gaussian (CSCG) random vector with mean x and covariance 
matrix S is denoted by CJ\f(x. X); and ~ stands for “distributed as”. C xxy denotes the space of x x y complex 
matrices. ||a;|| denotes the Euclidean norm of a complex vector x. For two real vectors x and y, x > y means 
that x is greater than or equal to y in a component-wise manner. 

II. System Model 

This paper considers the uplink communication in a multi-antenna C-RAN. As shown in Fig. |Tjc), the system 
consists of one BBU, N RRHs, denoted by the set Nf = {1, - - - ,N}, and K users, denoted by the set 1C = 
{1, ■ ■ • , A"}. It is assumed that each RRH is equipped with M > 1 antennas, while each user is equipped with one 
single antenna. It is further assumed that each RRH n is connected to the BBU via a digital error-free fronthaul 
link with a capacity of T n bits per second (bps). In the uplink, each RRH processes the signals received from all 
the users into digital bits and forwards them to the BBU via its fronthaul link. Then, the BBU jointly decodes the 
users’ messages based on the signals from all the RRHs. The details of our studied uplink multi-antenna C-RAN 
are given as follows. 

It is assumed that all the K users transmit over quasi-static flat-fading channels over a given bandwidth of B 
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RRH 1 



Fig. 2. Illustration of C-RAN with SCF applied at each multi-antenna RRH 

J The equivalent baseband complex symbol received at RRH n is then expressed as 

K 

y n = h n ,kVPk s k + z n = H n P*s + z n , n = 1, • • ■ ,N, (1) 

fc=l 

where s = [si, - - - ,sk] T with Sk ~ CJ\f(0, 1) denoting the transmit symbol of user k which is modelled as a 
CSCG random variable with zero-mean and unit-variance, P = diag(pi, • • • ,Pk) with pk denoting the transmit 
power of user k, H n = [h Ui i, • • • , h n) x] € C MxK with h n i- € C Mxl denoting the channel vector from user k 
to RRH n, and z n ~ CN{(). rr‘^1) denotes the additive white Gaussian noise (AWGN) at RRH n. It is assumed 
that s are independent over n. 

At each multi-antenna RRH, we consider the use of the SCF scheme to forward the received baseband symbol 

y n to the BBU via the finite-capacity digital fronthaul link. An illustration of the SCF scheme is shown in Fig. 

0 where each RRH first demodulates the signal received from each antenna to the baseband and processes the 

baseband signals from all antennas via spatial filtering, then conducts a scalar quantization over each output in 

1 It is worth noting that this paper studies the narrowband wireless channel, while the results obtained can be readily extended to broadband 
based C-RAN by viewing each frequency subchannel as one norrowband channel considered here and accordingly scaling down the fronthaul 
capacity T n by the number of frequency subchannels in use. 
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parallel, and finally forwards the total quantized bits to the BBU via the fronthaul link. Specifically, RRH n first 

implements spatial filtering to its received signal y n to obtain 

I< 

Un — ^nVn — ^n ^ ^ hiiyk^/PfoSk T" Vn^n — Vn,H n P^ S -(- Vn^m ^ — 1; ' ‘ ■ A", (2) 

k =1 

where V n = v„ \ , ■ ■■ . v n j^ r E C LxAI denotes the filtering matrix at RRH n with L < M denoting the 
(reduced) dimension of the output signal y n after filtering, which will be specified later. Then, a simple uniform 
scalar quantization is applied to each element of y n = [y n y, • • • , y n ,L] T at RRH n. Note that each complex symbol 
y n l can be represented by its in-phase (I) and quadrature (Q) parts as 

yn,l = y I n,l+3Vn,V h 0) 


where j 2 = — 1, and the I-branch symbol y J nl and Q-branch symbol y^ l are both real Gaussian random variables 
with zero mean and variance {Y^=iPk\^n ih n ,k\ 2 + <7 2 ||i; rij /|| 2 )/2. A typical method to implement the uniform 
quantization is via separate I/Q scalar quantization, the details of which can be found in lfT6l . After uniform scalar 

quantization, the baseband quantized symbol of y n is then given by 

K 

y n =y n + e n = VnYh n , k ^ k s k + V n z n + e n , n = l,---,N, (4) 

fc=i 


where e n = [e nj i, • • • , e n j,] r with e n j, 1 < l < L, denoting the quantization error for y n j with zero mean and 
variance q n j. Let D n j denote the number of bits that RRH n uses to quantize the I-branch or Q-branch of y n i, 

i.e., yl j or yz l . The quantization noise level q n j due to e n j for uniform quantization is then given by lfl6ll 

/ k \ 


Qn.l = 


Yl Pklv^ihnM 


k= 1 


2 + a 2 | 


V 


n, L 


>-2 D n 


if D n j. > 0, 


1 = L 


,L, n = 1, 


,N. (5) 


oo, if D n i = 0, 

Note that e n /s are independent over l due to independent scalar quantization for each element of y n , and also 

over n due to independent processing at different RRHs. As a result, the covariance matrix of e n is a function of 
P = [pi, ■ ■ ■ ,Pk] T , V n as well as D n = [D Ujl , • • • , D niL ] T , which is given by 


Qn(P, v n,D n ) = E[e n e%] = diag(g n ,i, • • • ,q n ,L)- 


( 6 ) 


Then, each RRH forwards the quantized bits to the BBU via the fronthaul link. The transmission rate in RRH 
n’s fronthaul link is expressed as llT6ll 

L 

T n = 2Bj2 D n,l, n = 1, - ■ ■ ,N. (7) 

i=i 


It is worth noting that if D n j = 0, then the symbol y n j at the spatial filter output dimension l is not quantized at 
RRH n and thus not forwarded to the BBU; as a result, q n j = oo according to (0). To summarize, our proposed SCF 



scheme is a two-stage compression method, where the output signal dimension at each RRH is first reduced from 
M to L (if M > L) by spatial filtering, and then further reduced by allocating the limited fronthaul quantization 
bits to only a selected subset of the L output dimensions with D n j > 0. 

The received signal at the BBU from all RRHs is expressed as 

K 

y = [yf,-- > Vn] T = V^h^^Sk + Vz + e, (8) 


k =1 

T . . . uT iT „ _ 


where V = diag(Vi,--- ,V N ), h (k) = [h lk ,-__ ,h Ntk ] 2 , z = [z[ , 

decode s k , we consider that a linear beamformingj is applied to y, i.e., 


and e = 


- r 


'1 5 


,e T N ] T . To 


h = w k y = w k Vh ik - ) y/p k s k + ^ j w k Vh^- )y /p j s j + w k Vz + w k e, k = l,---,K, (9) 

o+k 

where w k = [w kt i,--- , Wk,NL] T € C ArLxl is the receive beamformer for s k . Note that w k r n _^ L+ i = 0 if the 
Ith signal dimension at RRH n is not quantized, i.e., I) n j = 0. According to (O, the SINR for decoding s k is 
expressed as 


= _ Pk\w k Vh{k)? _ , i R 

E Pj\ w k Vh (j)\ 2 + (j2 \\ w k v \\ 2 + w^Q(p,{V n) D n })w k 

where Q(p,{V n ,D n }) = E[ee H } = diag(Q 1 (p, V ± , D ± ), ■ ■ ■ ,Q N (p,V N ,D N )). 


( 10 ) 


III. Problem Formulation 

In this paper, we aim to maximize the minimum SINR of all the users by optimizing the users’ power allocation, 
i.e., p, the spatial filter and quantization bits allocation at each RRH, i.e., V n and D n , Vn, and the receive 
beamforming vectors at the BBU, i.e., w k , \/k. Specifically, we aim to solve the following problem: 


Maximize 

7 ,p,{V n ,D n },{w k } 

Subject to 
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(id 


_ nWN) I 1 _ VJt 

E Pj\ w k vh (j)\ 2 + cr 2 \\w^v\\ 2 + wj*Q(p, {V n ,D n })w k 


( 12 ) 


Pk — Vfc, 


(13) 

L 



2 B ^ D n j < T n , 
1=1 

Vn, 

(14) 

D n ^i is an integer, 

Vn, l, 

(15) 


where 7 denotes the common SINR target for all the K users and P k denotes the transmit power constraint of 
user k. 

; It is worth noting that the subsequent results also hold similarly for the case with the non-linear successive interference cancellation 
scheme with any fixed decoding order, which is not considered in this paper due to its high complexity for implementation. 







Fig. 3. A diagram for illustrating the proposed algorithm for Problem GB¬ 
it can be observed that problem (fTTI) is a non-convex optimization problem since the design variables are 
complicatedly coupled in its constraint (fTTI) . Note that with given V n ’s and by setting I) n j = oo, Vra, l, problem 
(fTTI) reduces to the well-known SINR balancing problem via transmit power control and receive beamforming only. 
This problem has been efficiently solved by ||T91 . li20l based on the non-negative matrix theory lETfl . However, due 
to the new quantization noise term in the SINR expression, i.e., wj?Q(p, {V n . D n })wk, it can be shown that the 
algorithm proposed in fT9fl . If20fl cannot be directly applied to obtain the optimal power control and beamforming 
solution to problem (fTTI) even with given V n ’s and finite values of D n / s. As a result, with the additional variables 
Va s and integer variables D n ’ s, it is generally difficult to globally solve problem (fTTI) . 

In this paper, we propose an efficient algorithm to solve problem (fTTI) suboptimally, for which the main procedures 
are illustrated in Fig. [3] Specifically, in the channel training phase, all the users send orthogonal pilot signals to 
the RRHs with equal power P < mini<jt<x P}-. Let y p a 0 t. n € C ;Ux 1 denote the received signal at RRH n during 
the channel training phase. Then, each RRH n estimates the MIMO channel H n (assumed to be perfect) and 
computes the covariance matrix of its received signal (with the presumed knowledge of P and a 2 ), i.e., 

S n = E[y P i lot , n yp ilot}n \ = H n PH f + a 2 1 = PH n H* + a 2 1, n = 1, • • • , N. (16) 

Then each RRH independently determines its spatial filter V n based on S n . Next, each RRH n sends the effective 
channel V n H n = H n = [h nj i, • • • , h n ^] to the BBU. Then, based on its global CSI, i.e., H n ’s, n = 1, • • • , N, 
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the BBU computes p, D n ’s and w k s by solving the following reduced problem of problem (fill) : 


Maximize 

7 ,p,{D n },{w k } 

Subject to 


7 


( 17 ) 


Pk\w%h (k) \ 2 

- TT= - Tr - 77 - > 7, vk 

Pj\w%h(j )\ 2 + (J 2 \\w% V\\ 2 + wj^Q(p, {V n ,D n })w k 


(18) 


Pk < Pk, Vfc, 


(19) 

L 



2 B ^ D n i < T n , 

1 = 1 

Vn, 

(20) 

D n i is an integer, 

Vn, l, 

(21) 


where = [h lk ,--- ,h N ^ k ] T . Finally, the BBU sends D n to RRH n, n = 1,-- - ,N, via its fronthaul link 
for the implementation of scalar quantization, and p to one RRH (say RRH n), which then broadcasts p to all 
the users in the downlink. In the following section, we first show how each RRH computes its spatial filter in a 
distributed manner, and then discuss how the BBU solves problem dl 7 b . 

IV. Proposed Algorithm for Problem (fTTh 

In this section, we present details of the two-stage algorithm for solving problem dTTT > based on the diagram 
shown in Fig. [3] 


A. Distributed Spatial Filter Design at RRHs 

First, we show how to obtain the spatial filter at each RRH. It is observed from problem (fTTh that each user’s 
SINR at the BBU depends on V n ’s, n = 1, ■ ■ ■ . N, at all RRHs, which are also coupled with other design variables. 
As a result, it is difficult to obtain the jointly optimal spatial filters at all RRHs. Therefore, in the following we 
propose to design V n ’s in a distributed manner, where each RRH n computes its own spatial filter only based on 
its local information, i.e., S n in ([T6l ). First, each RRH n obtains the eigenvalue decomposition (EVD) of S n as 

S n = U^X n U n , n = !,-■■ ,N, (22) 

where A n = diag(A n q, • • • , A Hj m) € C MxM with A n q > ■ ■ ■ > A n ,M > 0 denoting the eigenvalues of S n , and 
U n = [u n> 1, • • • € C MxM is a unitary matrix with u n/ni denoting the eigenvector of S n associated with 

the eigenvalue A n> m, tn = 1, • • • , M. Notice that U n can also be obtained as the singular vectors of H n from 
its singular value decomposition (SVD). Note that since the rank of S n is rank(^ n ) = min (M,K), A n ^ m = 0 if 
min(M, K ) < m < M, and thus each RRH only needs to extract the first L = min(M, K ) user signal subspace 
by setting its spatial filter as 


Vn — Un — [^n,l j > ^n,L] j tl — !,•••, N. 


( 23 ) 
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As a result, the equivalent channel from all users to RRH n after filtering is given by 

H n = U n H n , n = 1, • • • ,N, (24) 

which is sent to the BBU in channel training as discussed in Section Hill 

The reasons to set the spatial filter at each RRH as in ( 123T ) are as follows. First, if M > K, U n denoises 
the received signal before quantization by extracting the user signal subspace from the received signal plus noise 
without loss of information lf22l . lf23l . Second, if M > K, the signal dimension is reduced from M to L = I\ at 
each RRH n, which simplifies the subsequent uniform scalar quantization at each dimension. Last but not least, 
with this spatial tiler applied at RRH n in data transmission, the filter output y n given in d2]) has the following 
covariance matrix, 

E[y n y% } = U n H n PH* + a 2 1, (25) 

which is (approximately) diagonal if the transmit power of users, i.e., pk s, are (closely) equal, i.e., the elements 
in y n are approximately uncorrelated. As a result, independent scalar quantization at each output dimension of y n 
as shown in © becomes more effective. 


B. Joint Optimization at BBU 

Next, given the spatial filter at each RRH V n = U n , Vn, as in (l23l) . we jointly optimize users’ power allocation, 
RRHs’ quantization bits allocation, and BBU’s receive beamforming by solving problem (fT71) . 

We propose to apply the alternating optimization technique to solve problem (fTTT ). Specifically, first we fix the 
quantization bits allocation D n i = D n j, Vn, l, in problem ([T71 ) and optimize the users’ power allocation and BBU’s 
receive beamforming solution by solving the following problem. 


Maximize 

7 ,p,{w k } 


Subject to 


7 

_ Pk\w^h {k) \ 2 _ 

Yj Pj\ w k h(j )\ 2 + cr 2 ||m fc || 2 +w%Q(p,{U n ,D n })w k ~ 7 ’ 

j¥=k 

Pk < Pk, Vfc. 


Vfc, 


(26) 

(27) 

(28) 


Let p = [p \, • • ■ ,pk] T and m/.’s denote the solution to problem (l26l ) . Then, we fix the users’ power allocation and 
BBU’s receive beamforming solution as p = p, = ijbk s, respectively, and optimize the RRHs’ quantization 
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bits allocation by solving the following problem. 


Maximize 

7 ,{D„} 

Subject to 


7 

_ Pk\w k h(k )\ 2 _ 

E + °' 2 ll™fcll 2 + w k Q{p,{U n ,D n })w k ~ 7 ’ 

j¥=k 

L 

2 B'Yh D n,l<T n , Vn, 

i=i 


D n i is an integer, Vn, l. 


VA:, 


(29) 

(30) 

(31) 

(32) 


Let Dn s denote the solution to problem ( l29l ). The above procedure is iterated until convergence. In the following, 
we show how to solve problems (l26l) and (l29l ). respectively. 

1. Solution to Problem ( 1261 ) 

First, we solve problem (1261 ). As mentioned in Section [Dll due to the additional term accociated with the 
quantization noise, i.e., w k Q(p,{U n ,D n })w k s, in the SINR expression, which is coupled with p and w k s, 
the algorithm proposed in |T9l and EOl cannot be applied to obtain the optimal solution to problem ( l26l ). In the 
following, we propose to apply the bisection method |[24l together with the celebrated fixed-point method lfl 8 l 
to globally solve problem (l26l) . Specifically, given any common SINR target 7 for all the users, we solve the 
following feasibility problem. 


Find p,{w k } 
Subject to 


p k \w^h (k) \ 


E Pj\ w k^(j)\ 2 + cr 2 ll™fcll 2 + w%Q(p,{U n ,D n })w k 


> 7 , VA:, 


(33) 

(34) 


Pk < Pk, VA:. (35) 

The feasibility problem (l33l) is still a non-convex problem. In the following, we show that the feasibility of 
problem (l33l) can be efficiently checked by solving a sum-power minimization problem without the individual 


power constraints in (I35I) . The sum-power minimization problem is formulated as 

I< 

Minimize > p k 

P '{ w h 


Pk\w%h (k) | 


Subject to 


> 7 , VA:. 


E Pj\ w k h (j)\ 2 + cr 2 ll w fcl | 2 + w^Q(p,{U n ,D n })w k 

j^k 

Note that given any power allocation p, the optimal linear receive beamforming solution to problem 


minimum-mean-square-error (MMSE) based receiver given by 

H 


-1 


w 


MMSE 


+cr 2l + Q(p, {U n ,D n }) ] h (k) , A: = 1 , - - - ,K. 

j^k 


(36) 

(37) 
is the 

(38) 
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With the above MMSE receivers, the SINR of user k given in (flOl ) reduces to 

-l 


Ik = Pkt\ k) | + a 2 1 + Q(p, {U n , D n }) ] h (k) , k = 1, ■ ■ ■ , K. 


-H 


As a result, problem (1361) reduces to the following power control problem. 



K 

Minimize 

P 


k= 1 

Subject to 

Pkf\ k) | 

) = [ht P)," 

• ,Ik(p)] 


■H 


xl 


Ik(p) = 


7 


k = !,■■■ , K. 


■H 


■ H 


h {k) E Pj h (j) h (f) + (j21 + Q(P. i U n, Dn}) h {k) 
\3^k I 


Then, problem (1401 reduces to the following problem. 


K 


Minimize p k 

P k =l 


Subject to p > ^(P)- 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


Lemma 4.1: I(p) given in (l42l) is a standard interference function f! 8 l . In other words, it satisfies 1. I(/p) > 0 
if p > 0; 2. if p > then J(p) > I(p')\ and 3. Va > 1, it follows that al{p ) > /(op). 

Proof: Please refer to Appendix |A| ■ 

Based on Lemma 14.11 we have the following corollaries. 

Corollary 4.1: If 7 is not achievable for all the K users in problem (l40l) . the iterative fixed-point method 
pL+i) — J(pW) with the initial point p(°) = 0 , where pW denotes the power allocation solution obtained in the 
ith iteration of the above fixed-point method, will converge to p oc. 

Proof: Please refer to Appendix |B] ■ 

Corollary 4.2: If 7 is achievable for all the K users in problem (l40l) . the iterative fixed-point method p^ l+l ) = 
J(p(*)) will converge to the optimal solution to problem (|43T ) with any initial point p (<>> > 0 . 

Proof: Please refer to llT8l Theorem 2], ■ 

Corollary 4.3: If 7 is achievable for all the K users in problem (l40l) . the optimal power solution p* to problem 
(1431 ) is component-wise minimum in the sense that any other feasible power solution p' that satisfies (1441 ) must 
satisfy p' > p*. 


Proof: Please refer to Appendix 0 
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According to Corollaries I4.1H4.31 the feasibility of problem (1331) can be efficiently checked as follows. First, we 
apply the fixed-point method in Corollary 14.11 to solve the sum-power minimization problem (l43l) with the initial 
point p <ir> = 0. If the obtained solution, denoted by p*, is infinity (unbounded), then 7 cannot be achieved by all 
the K users simultaneously even without the individual power constraints. Therefore, problem (l33l) is not feasible. 
Otherwise, if p* is of finite value, 7 can be achieved by all the users without the individual power constraints. In 
this case, we check whether p* satisfies all the given individual power constraints. If p* does not satisfy all the 
individual power constraints, according to Corollary 14.31 all the power solutions that satisfy the SINR constraints 
cannot satisfy the given individual power constraints. As a result, problem (|33T ) is not feasible. Otherwise, if p* 
satisfies all the individual power constraints, it is a feasible solution to problem (f33T >. i.e., problem ( |33T ) is feasible. 

Let 7* denote the optimal value of problem (l26l) . The algorithm to solve problem ( 1261 ) based on the bisection 
method is summarized in Table |U Note that in Step 1, 7 max can be set as 7 max = max] <^.<77 77, where 77 is 
obtained by setting />/,. = F\ and 77 = 0 , Vy / k, in (l39l) . 

TABLE I 

Algorithm!) Algorithm for Problem d26l > 


1. Initialize 7 min = 0, 7 max > 7*5 

2. Repeat 

a) Set 7 = 2 max± 7 min. 

b) Solve problem 1431 by the fixed-point method as in Corollary 14.11 with initial point = 0; 

c) If the obtained solution p * is infinity or does not satisfy all the given individual power constraints in 1281 . set 7 max = 7 ; 
otherwise, set 7 m i n = 7 ; 

3. Until 7 max — 7 min < e, where t > 0 is a small constraint to control the algorithm accuracy; 

4. Calculate loLs by substituting p* into ( 1381 . 


2. Solution to Problem ( 1291 ) 

Next, we consider problem (1291 ). Note that the quantization noise power given in ©, i.e., q n j, is a continuous 
function over quantization bits allocation D n [ (assuming it is continuous), except when I)„ i = 0. To deal with 
this issue, we approximate the quantization noise power with the following continuous function of I) n .i, i.e., 

q n ,i = 3 (^2Pk\Vn,iK,k\ 2 + °' 2 |l w n,/|| 2 ^ 2 -2A *' ! , D n q> 0, Vra, Z. (45) 

The difference of (l45l) from © is that when D n j = 0, q n j equals Pk\ v n An,k\ 2 + cr 2 ||i; n ^|| 2 ) instead of 

infinity. This indicates that even with I) ri j = 0, some quantized information at the Zth signal dimension of RRH 
n is forwarded to the BBU. However, since 3(^(■ =r 1 Pt\ v i^-n.k\ 2 + fjl \\ v n,i\\ 1 ) is generally much larger than the 
signal power of any user k, i.e., Pk\v„ ih n! k\ 2 ’s, its effect on the user SINR given in (fTOl ) is negligible. Therefore, 
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we use (l45l) to approximate the quantization noise power for solving problem (l29l ) . Note that when the actual SINR 
is computed using (fTOl) . the quantization noise power is set to be infinity if the obtained D n i is zero. 

With the approximation (l45l ). problem (l29l ) is still challenging to solve due to the integer constraints for D n / s. 
Let problem (1291-NoInf) denote the relaxation of problem (l29l ) without the integer constraints in (l32l) . In the 
following, we first solve problem (1291-NoInt). the solution to which may not satisfy all the integer constraints. 
Then, we propose an efficient algorithm to obtain a set of integer solutions for all D n /s based on the solution of 
relaxed problem (1291-NoInf). 

Similar to problem (l26l ). problem (1291-NoInt) can be globally solved by the bisection method. Given any common 
SINR target 7 for all the users, we need to solve the following feasibility problem. 


Find {£>„} 
Subject to 


p k \w%h {k) \‘ 


Y Pj\ w k h U)\ + 0 " ll^fcll + w^Q{p,{U n ,D n })w k 
j¥=k 

L 

2 B^D n j<T n , Vn. 

z=i 


> 7 , VAs, 


Note that we have 


N L 


N L 


w k Q{Pi {Un, D n })w k = ^2^2q n j\w k . 


(n— l)L+m| 


y: yi Qn,i,k 2 


—2D„ 


n=l 1=1 


n= 1 1=1 


where Wi , denotes the jth element of w, , 1 < i < K, 1 < j < NL, and 


i,l,k ^\^k,(n—l)L+l\ ^ Pil^ri,i\ T Z. 


VAs, 


(46) 

(47) 

(48) 


(49) 


(50) 


Note that 9 n i k can be interpreted as the effective quantization noise power due to the Zth quantized dimension at 
RRH n in decoding s k at the BBU. Then, problem (1461 is equivalent to the following problem. 


Find {D n } (51) 

N L - | - Hu 1 2 

Subject to ^T^2d n ,i } k2~ 2Dn ’ 1 < —— LJY. - ^p^w^h^l 2 - a 2 \\w k \\ 2 , Vfc, (52) 

n= 1 1=1 ' j^k 

L 

2 Bj2 D n,i<T n , Vn. (53) 

l=i 

Problem (15Tb can be shown to be a convex feasibility problem, and thus can be efficiently solved via the interior- 
point method l24l . Let 7 * denote the optimal value of problem (1291-NoInt). The algorithm for problem (1291-NoInt) 
is then summarized in Table 1TT1 Note that in Step 1, 7 max can be similarly set as in Algorithm Q] 

It is worth noting that the solution obtained by Algorithm |TTj may not satisfy all the integer constraints given by 
(l32l) in problem (l29l ). In the following, we show how to obtain a set of integer solutions {D n } for problem (l29l ) 
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TABLE II 

Algorithm|TT1 Algorithm for Problem (|29INoInt) 


1. Initialize 7 min = 0, 7 max > 7*5 

2. Repeat 

a) Set 7 = Tmax+7min . 

b) Solve problem ( 151b by the interior-point method. If problem ( 151b is feasible, set 7 m in = 7 ; otherwise, set 7 max = 7 ; 

3. Until 7 m ax — 7min < e, where e > 0 is a small constraint to control the algorithm accuracy. 


based on the solution of relaxed problem (1291-NoInt). Similar to lfl6ll . we propose to round each I) n j to its nearby 
integer as follows. 

A f [D n i\, if D n i — |_-Dn,zJ < a n , 

D n ,i = < n = I,-- ,N, l = l,--- ,L, (54) 

[ | D n j |, otherwise, 

where 0 < a n < 1, Vn, and [x \, [:/;] denote the maximum integer that is no larger than x and the minimum 
integer that is no smaller than x, respectively. Note that we can always find a feasible solution of D n f s by simply 
setting a n = 1, Vn, in (l54l) since in this case we have YliL 1 D n ,i < Ef =1 D n , < T n /2B. In the following, we 
show how to find a generally better feasible solution by optimizing a n ’ s. Since D n /s decrease with a n ’s, as a:„’s 
become smaller, the resulting D n f s from (l54l) achieve higher SINRs for all the users, but the fronthaul constraints 
for the RRHs become more difficult to be satisfied. Hence, we propose to apply a simple bisection method to find 
the optimal values of a n ’ s, and then substitute a„’s into (l54l) to obtain D n /s. The algorithm is summarized in 
Table [HD 

TABLE III 

Algorithm IIIII Algorithm to Find Feasible Integer Solution D„ t i ’s for problem 


1. Initialize ot™ m = 0, a“ ax = 1, Vn; 

2. For n = 1 : N, repeat 

„ min 1 „ max 
o . oc„ +Q!„ 

cl. o6t OLxi — 2 ’ 

b. Substitute a n into 154b . If D n / s, VI, satisfy 2 B N.fl] D nt i < T n , set a™ ax = a n \ otherwise, set Q™ ln = a n \ 

3. Until a™ ax — a™ ln < e, Vn, where e > 0 is a small constraint to control the accuracy of the algorithm. 


3. Overall Algorithm for Problem rtf 71 ) 

After problems (l26l) and ( 129b are solved by Algorithm ID and Algorithms [TT] and |TID respectively, we are ready to 
propose the overall algorithm for problem (fTTT) with given V n = U n , Vn, based on the alternating optimization 
technique, which is summarized in Table |IV] The algorithm starts with equal quantization bits allocation for all 
the spatial filter output dimensions at each RRH, i.e., I) n j = [T n /(2BL)\, Vn, l. Note that the proposed algorithm 
terminates in either of the following two cases: 1. 0 < 7 W — 7 V- 1 ) < 7 i.e., the minimum SINR of users cannot 
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be improved above the positive threshold; and 2. 7 W — 7 b -1 ) < 0, i.e., the minimum SINR obtained at the current 
iteration is even reduced compared with that at the previous iteration. The second case may occur since problem 
(l29l ) is generally not globally solved by Algorithms HT1 and Hill due to the integer constraints in (l32l) . 

TABLE IV 

Algorithm[Tv1 Algorithm for Problem (1171) with given V n = U„ 


1. Initialize: Set D^\ = \T n /(2BL) J, Vn, l, = 0, and i = 0; 

2. Repeat 

a. i = i + 1 ; 

b. Obtain p b) and [wf 1 } by solving problem ( 1261 ) with D n y = Vn, Z, using Algorithm Q] 

c. Obtain {Dr , 1 } by solving problem 1291 ) with p = p hi and Wk = w ^\ Vfc, using Algorithms [III and Hill 

3. Until 7 b) — w here 7 b) denotes the objective value of problem 1171 achieved by p^\ {wf 1 } and {Df }, and e > 0 is 

a small constraint to control the accuracy of the algorithm. 


V. Numerical Results 

In this section, we provide numerical results to verify our results. In the following numerical examples, the 
bandwidth of the wireless link is B = 10MHz, while the path loss model of the wireless channel is given as 
30.6 + 36.7log 10 (72) dB, where d (in meter) denotes the distance between the user and the RRH. The transmit 
power constraint for each user is Pk = 23dBm, Vfc. The power spectral density of the background noise at each 
RRH is assumed to be —169dBm/Hz, and the noise figure due to the receiver processing is 7dB. Furthermore, it 
is assumed that all the RRHs possess the identical fronthaul capacity, i.e., T n = T, Vn. 

A. Effectiveness of Distributed Spatial Filtering at RRHs 

First, we illustrate the effectiveness of the proposed RRHs’ distributed spatial filtering via SCF. In this example, 
there are iV = 4 RRHs and I\ = 8 users randomly distributed in a circle area of radius 500m. Besides the 
proposed spatial filter design given in (l23l) at each RRH, we also consider the following schemes for performance 
benchmark: 1. Matched filtering, i.e., V n = H^, Vn; 2. Zero-forcing (ZF) filtering, i.e., V n = 

Vn; and 3. Without spatial filtering, i.e., V n = I, Vn. Note that for the matched and ZF filtering, like the proposed 
design, the signal dimension is reduced from M to K at each RRH, since in this example we set K < M. 

Fig. 0] shows the maximized minimum (max-min) SINR of all the users achieved by Algorithm HVl with different 
spatial filtering designs at the RRHs, where each RRH is equipped with M = 10 antennas and the per-RRH 
fronthaul capacity T varies from T = O.IGbps to T = 1.5Gbps. It is observed that the max-min SINR achieved 
with the proposed spatial filtering in (1231) based on the EVD of the covariance matrix S n (or SVD of H n ) is higher 
than that achieved by the other benchmark schemes for all values of fronthaul capacities. This observation also 
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Fig. 4. Performance comparison of different distributed spatial filtering designs at RRH versus per-RRH fronthaul capacity. 

holds in other simulation setups, which are omitted due to the space limitation. The performance gain is explained 
as follows. For matched filtering, the dimension of the spatial filter output y n given in © is reduced from M 
to L = K, but its elements are still correlated at each RRH n, which makes the subsequent scalar quantization 
ineffective. On the other hand, with ZF filtering, although the signal dimension is also reduced and the elements of 
y n are decorrelated in general, the noise effect becomes more severe after filtering and quantization, which degrades 
the joint decoding SINR at the BBU. Furthermore, with V n = I, Vn, our proposed SCF scheme reduces to the 
antenna selection based quantization where the received signals from some antennas are not forwarded to the BBU 
if the quantization bits allocation obtained from Algorithm [IV] yields D n , m = 0 for antenna m at RRH n. With 
limited per-RRH fonthaul capacity, this scheme is not effective as substantial user signal information is discarded 
at each RRH due to antenna selection, even with optimized quantization bits allocation. Last, it is observed that as 
the per-RRH fronthaul capacity increases to 1.5Gbps, the max-min SINRs achieved with all spatial filter solutions 
converge to the same value. This is because when the fronthaul capacity is large enough, the quantization error 
becomes negligible as compared to the receiver noise at each RRH. 

Fig. [5] shows the max-min SINR of all the users achieved by Algorithm [TV] with different spatial filtering 
designs at the RRHs versus the number of antennas at each RRH, where the per-RRH fronthaul capacity is fixed 
as T = 0.5Gbps. It is observed that the proposed filter outperforms both matched filter and ZF filter for all values 
of M. Moreover, it is also observed that as M increases, the performance gap between the schemes with versus 
without filtering is enlarged. This is because when M is large, we need to compress the signals received at each 
RRH more significantly by spatial filtering before scalar quantization, i.e., with L <C M. 
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Fig. 5. Performance comparison of different spatial filter design at RRH versus per-RRH number of antennas. 

B. Performance Gain of Joint Wireless and Fronthaul Resource Allocation 

Next, we show the performance gain of joint optimization for wireless and fronthaul resource allocation. In 
addition to the proposed Algorithm |WJ we consider the following three benchmark schemes for performance 
comparison. 

• Benchmark Scheme 1: Optimized quantization bits allocation without power control. In this scheme, all 
the users transmit at their maximum power, i.e., Pk = Pk, Vfc, while the BBU iteratively updates its MMSE 
based receive beamforming as given in (1381) and RRHs’ quantization bits allocation via Algorithms UTI and Hill 

• Benchmark Scheme 2: Optimized power control with equal quantization bits allocation. In this scheme, 
each RRH equally allocates its fronthaul capacity to all the L dimensions of the output signal after spatial 
filtering, i.e., D n j = T n /{2BL), \/l,n. Then, the BBU computes the users’ power allocation p as well as its 
beamforming solution wfs using Algorithm Jl 

• Benchmark Scheme 3: Equal quantization bits allocation without power control. In this scheme, all the 
users transmit at their maximum power, and each RRH equally allocates its fronthaul capacity to all the L 
dimensions of the filter output signal, while the BBU computes the resulting MMSE receive beamforming. 

Fig. [6] shows the performance comparison of different schemes versus per-RRH fronthaul capacity T under the 
same setup as for Fig. [4] It is observed that the proposed joint optimization of wireless and fronthaul resource 
allocation, i.e., AlgorithmII VI achieves higher max-min SINR over Benchmark Schemes 1-3 with separate resource 
allocation. Moreover, it is observed that under small and moderate per-RRH fronthaul capacity, Benchmark Scheme 
1 performs closer to Algorithm [[V] as compared to Benchmark Schemes 2 and 3, which equally quantize all the 
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Fig. 6. Performance comparison of different resource allocation schemes in the case of small number of users. 

TABLE V 

Number of Quantized Dimensions at Each RRH 


RRH Index 

1 

2 

3 

4 

No. of Quantized Dimensions (T=0.lGbps) 

3 

3 

3 

3 

No. of Quantized Dimensions (T=0.5Gbps) 

6 

7 

5 

7 

No. of Quantized Dimensions (T=1.0Gbps) 

8 

8 

8 

8 


signal dimensions at each RRH. In other words, in this example, most of the joint resource allocation gain is due to 
joint quantization bits allocation at RRHs. Note that in this example, the total number of antennas from all RRHs 
is MN = 40, which is much larger than that of the users K = 8. Since there are sufficient spatial dimensions to 
separate the users’ signals, power control becomes less effective than quantization bits allocation. Table lYl shows 
the number of quantized dimensions with I) n j > 0, / =1. • • • , L, at each RRH n when the per-RRH fronthaul 
capacity is T = 0.1, 0.5, lGbps, respectively. It is observed that when T = O.IGbps, each RRH should quantize 
3 out of L = min(M, K) = 8 signal dimensions. As the fronthaul capacity increases, more signal dimensions 
can be quantized. Particularly, when the fronthaul capacity is sufficiently large, i.e., T = lGbps, all the 8 signal 
dimensions after spatial filtering should be quantized at all RRHs. 

Fig. |7] shows the performance comparison for the case of large number of users with I\ = 20. It is observed that 
Algorithm [IV] still outperforms Benchmark Schemes 1-3. However, different from the case with smaller number 
of users (K = 8) shown in Fig. [6] it is observed that with K = 20, Benchmark Scheme 2 performs better than 
Benchmark Schemes 1 and 3. This is because with more users but fixed total number of spatial degrees of freedom 
MN = 40, the inter-user interference is more severe and thus power control becomes more effective. 
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Fig. 7. Performance comparison of different resource allocation schemes in the case of large number of users. 

C. Multi-Antenna C-RAN versus Massive MIMO 

As discussed in Section |]J besides multi-antenna C-RAN considered in this paper, massive MIMO and single¬ 
antenna C-RAN arc two promising techniques proposed for 5G wireless networks. Therefore, an interesting as 
well as important question we seek to address in this subsection is as follows: given a total amount of antennas 
to be deployed, should we equip them in one single BS, i.e., massive MIMO, or distribute them over a given area 
by connecting to the BBU via finite-rate fronthaul links, i.e., C-RAN? Moreover, if C-RAN is preferred, what is the 
optimal antenna deployment solution, i.e., single-antenna C-RAN versus multi-antenna C-RAN? Intuitively, if more 
single-antenna RRHs are deployed in the network, with higher probability each user can be served by one or more 
nearby RRHs with strong channel conditions. However, with multi-antenna RRHs, we can efficiently perform SCF 
at each RRH to better utilize the limited fronthaul capacity given for each RRH. In the following, we provide a 
case of study to show the advantage of multi-antenna C-RAN over its two extreme counterparts: massive MIMO 
with all the antennas deployed at one BS and single-antenna C-RAN with only one antenna at each RRH, i.e., 
M = 1. 

To make a fair comparison, we assume that there are in total M antennas to serve K users in a given area. 
Specifically, for the massive MIMO system, we assume that there is only one BS which is equipped with all M 
antennas, while for the C-RAN, we assume that there are N, 1 < N < M, RRHs each equipped with M/N 
antennas^ Moreover, we assume there is a total fronthaul capacity for all RRHs in C-RAN, denoted by T, which 
is a constant regardless of the number of RRHs N. Note that this is a valid assumption since in practice the signals 

’In the case that M/N is not an integer, we can simply let some RRHs possess \_M/N\ antennas, while the others possess \M/N~\ 
antennas such that the total number of antennas is equal to M. 
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Fig. 8. Performance comparison between multi-antenna C-RAN and massive MIMO. 

from different RRHs in the same area are usually first multiplexed at a local hub, which then forwards the signals 
to the BBU via a fronthaul link with capacity T. For simplicity, we set T n = T / N for all the N RRHs in the case 
of C-RAN. 

Fig. [8] shows the max-min SINR performance comparison between massive MIMO versus C-RAN with different 
number of RRHs and fronthaul sum-capacity. The total number of antennas is M = 50. Moreover, there are 
K = 20 users randomly located in a circle area of radius 700m. For the massive MIMO, the BS is located in the 
center of the circle, while for the C-RAN, the RRHs are randomly located in the circle. Note that since the BS in 
massive MIMO system directly decodes the user messages, the optimal power allocation for users and decoding 
beamforming vectors at the BS can be obtained by Algorithm Q] with D n i = oo, Vn, l. It is observed from Fig. [8] 
that if the fronthaul sum-capacity is sufficiently large, i.e., T = 3Gbps or 5Gbps, the max-min SINR in C-RAN 
is larger than that in massive MIMO when the number of RRHs N > 4. This is because the densification gain 
due to larger number of RRHs dominates the performance of C-RAN. However, if the fronthaul sum-capacity is 
limited, i.e., T = lGbps, the performance of C-RAN with any N is inferior to that of massive MIMO, since 
quantization errors dominate its performance. Moreover, for the case of single-antenna C-RAN with N = 50 
single-antenna RRHs (which is not shown in Fig. [8] for brevity), the achieved max-min SINRs for the cases of 
T = 1,3,5Gbps are 2.25, 10.23 and 22.08, respectively, which are significantly lower than the maximum SINRs 
achieved in multi-antenna C-RAN with the corresponding optimal number of RRHs shown in Fig. [8] This is 
because in the case of single-antenna C-RAN, the received signals at nearby RRHs are highly correlated, and as 
a result independent scalar quantization at these RRHs without spatial filtering is inefficient in utilization of the 
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given fronthaul sum-capacity. Last, it is observed that the achieved max-min SINR in multi-antenna C-RAN first 
increases with the number of RRHs, but then decreases after a certain optimal number of RRHs for all three cases 
of T = 1,3,5Gbps. This validates the effectiveness of multi-antenna C-RAN by more flexibly balancing between 
the performance and fronthaul trade-off. As T increases, it is observed that more RRHs should be deployed in 
C-RAN to exploit better channels from all users. 

VI. Conclusions 

This paper considers a flexible antenna deployment design for C-RAN termed multi-antenna C-RAN and proposes 
a new “spatial-compression-and-forward (SCF)” scheme for efficient and low-complexity processing at each RRH 
in the uplink multiuser communication. With the proposed distributed spatial filters at RRHs, a joint optimization 
across the wireless transmission, the fronthaul quantization and the decoding at the BBU is performed to maximize 
the minimum SINR of all the users. Our results show that the proposed SCF scheme with joint resource allocation 
achieves significant performance gains over the conventional “quantize-and-forward” based single-antenna C-RAN 
as well as massive MIMO, thanks to the flexibly optimized antenna distribution in multi-antenna C-RAN. 


Appendix 


A. Proof of Lemma \4.1\ 


-H 


First, since Pjh(j)h(j) +& 2 I + Q(p,{U n ,D n }) >- 0, fk, we have I(p) > 0 if p > 0. The first 

_ / 

property in Lemma |4~TI is proved. 

Next, if p > p', then we have 


^2pjh(j)h%j) + (T 2 1 + Q(p , {U n , D n }) 




+^ I + Q(p',{Un,D n }), Vfc. 

Lemma A.l: [21, Corollary 7.7.4] If A y- 0, B y 0, and A — B V 0, then A -1 — B 1 A 0. 
According to Lemma IA. 1 1 it follows that 

■H 


^2Pjh(j)h {j) + a 21 + Q{p,{U n ,D n }) 


-l 


=< I Y,p'Aj) h u) + (j2/ + Q(p'. {Un, D n }) 

\j¥* 

As a result, we have I (p) > I(p') if p > p'. The second property in Lemma 14.11 is thus proved. 


, VL 


(55) 


( 56 ) 
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Last, if a > 1, then we have 


+(r 2 I + Q{p,{U n ,D n }) 


yjMk 


y ^2 a Pjhti)h(j) + o 2 I + Q{ap, {U n , D n }), Mk. 
j¥=k 

According to Lemma |4~T1 we have 

— j 'y -f- o I + Q(p, {U n , D n }) 




-1 


A | ^2 a Pjh(j)h(j) + a2j + Q( a P, { U n, Dn}) I , V/:. 


As a result, it follows that Va > 1, al{p) > I(ap). The third property of Lemma 14. 1 1 is thus proved. 
To summarize, I(p) given in (l42l) is a standard interference function. Lemma |4~T1 is thus proved. 


(57) 


(58) 


B. Proof of Corollary 14.21 

If the initial point of the fixed-point method is p^°) = 0, then we have p ^ = J(p(°)) > 0 = p(°\ Moreover, 
according to the second property of Lemma |4~T1 if pW > p^~ l \ then p(* +1 ) = 7 (pW) > I( y p^ 1 '>) = pW. As a 
result, with the initial point p = 0, the resulted power solution increases after each iteration, i.e., p^ < p ^ < 
p( 2 ) < • • •. Suppose that at the ith iteration, the above procedure converges to a finite power solution pW. Then 
we have ph +1 ^ = pW = J(pM). As a result, pW can satisfy the SINR constraints (1441 ) . which contradicts to the 
fact that 7 cannot be achieved by all the users. To summarize, the power solution by the fixed-point method with 
P'°) = 0 will increase to an infinity power allocation. Corollary 14.21 is thus proved. 


C. Proof of Corollary 14.71 

Let p(°) denote a feasible solution to problem (l43k Since it satisfies (l44l ). we have p ^ > I(p (°)) = pW. 
In the following, we show that pi l+l ) < pW holds Mi > 1 by induction. Suppose pM < p(* _1 (. Then we have 

(a) 

pfi+i) _ < J(pfi -1 )) = pW, where (a) is due to the second property in Lemma |4~T1 Since we have 

pi 1 ) < p(°) ) we thus have pi l+1 ) < pW holds Vi > 1. In other words, given any feasible power solution p’ as the 
initial point, a monotonic convergence can be guaranteed by the fixed-point method, i.e., p’ > p*. Corollary 14.31 
is thus proved. 
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